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I discuss low-energy lines in gamma-ray bursts. The process of de- 
convolving gamma-ray spectral data and the steps needed to demon- 
strate the existence of a line are explained. Previous observations and 
the current status of the analysis of the BATSE data are described. 



^ Spectral lines are highly informative of the conditions of the region in which 
they are created. Unfortunately, the observational and theoretical under- 
^ standing of lines in gamma-ray bursts is confused at best. Observations made 
^ ^ ■ with Ginga provided strong evidence for absorption lines, yet BATSE has not 
Q \ confirmed their existence. Confirmation of the existence of lines and a bet- 
■ ter characterization of their properties might provide the essential clue to the 
mystery of gamma-ray bursts. 

In this paper I will not tell you whether lines exist, instead I will explain 
what analysis needs to be done to demonstrate the existence of a line in a 
gamma-ray spectrum to assist you in judging the evidence for yourself. I will 
K^" \ briefly review selected observations and the status of the analysis of BATSE 
data. I will concentrate on observational issues regarding spectral features 
below 100 keV. The views expressed are my own. 



OO 



\Q DETECTORS AND GAMMA-RAY INTERACTIONS 

ON ■ 

i The detection of gamma-ray lines is much more difficult than the detection 
Oh of optical lines, primarily because there is no one-to-one relationship between 
Q . the energy of an incident gamma-ray photon and the energy measured by the 
I detector, which is called the "energy loss" . Because of the difference between 
. the energy of the incident photon and the energy loss, detected events are 
referred to as "counts" . Secondary problems are the poor signal-to- noise and 
signal-to-background ratios prevalent in gamma-ray astronomy. 

The examples in this paper are from data collected with the Spectroscopy 
Detectors (SDs) of BATSE, but the concepts are true for all gamma- 
' ray detectors — regardless of their energy resolution, the photon interaction 
physics is similar. The SDs are 12.7 cm diameter by 7.6 cm thick crystals of 
Nal(Tl) scintillator, each viewed by a photomultiplier tube (PMT) of the same 
diameter. When a gamma-ray interacts in the Nal crystal, a fraction of the re- 
sulting ionization energy is converted into scintillation light and measured by 
the PMT. In the energy range of interest, 10 keV to a few MeV, gamma-rays 
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FIG. 1. The Nal gamma-ray mass attenu- 
ation coefficient fi: soUd fine: total; dotted: 
photoelectric; dashed: Compton; dot-dash: 
pair production. The fractional transmission 
is exp{—fix), where x is the quantity of Nal 
traversed in g cm"'^. The data are from ref. 



interact with detector matter primarily by three processes ( Fig. |l|) : photo- 
electric absorption, Compton scattering, and pair production ( p"4|Jl^j20| ) . 

In Nal, the photoelectric process dominates below « 200 keV. In this process 
the gamma-ray is completely absorbed by an atomic electron. The interaction 
is most likely to occur with the inner-most electron that the gamma-ray has 
sufficient energy to ionize. Usually, an atomic cascade yields fluorescent X- 
rays, which will be photoelectrically absorbed unless they escape the crystal. 
Sometimes an Auger electron is ejected instead. The left curve of Fig. || illus- 
trates a simple case, the detected energy losses expected from 20 keV photons. 
The fluorescent X-rays have a higher cross-section because of their lower en- 
ergy (Fig. P, so that they are unlikely to escape and thus the incident energy 
will be totally absorbed, leading to a simple observed spectrum. Broadening 
of the spectral resolution occurs, primarily due to Poisson fluctuations in the 
number of photoelectrons produced in the PMT ( |20| ) . 

A 50 keV photon is also most likely to interact via the photoelectric process, 
but the resulting count spectrum is more complicated (Fig. ^ bold curve). 
The inner-most shell of an iodine atom accessible to the 20 keV photon of the 
previous example is the L-shell. The K-shell becomes accessible to photons 
with energies above the K-shell binding energy of 33.17 keV, resulting in a 
cross-section increase (Fig. |^) . A 50 keV photon will typically ionize a K-shell 
electron, resulting in a L, M or N to K shell transition and an X-ray with an 
energy between 28.3 and 33.0 keV (|2^). Sometimes this photon will interact 
by the photoelectric process and the total energy of the initial photon will 
be collected, resulting in the full-energy absorption peak at the right in the 
bold curve of Fig. |[ However, the fluorescence X-ray is below the K-edge 
and has a lower cross-section (Fig. |l|) and thus a high probability of escaping 
the crystal, resulting in incomplete absorption of the energy of the incident 
gamma-ray (the left peak of the bold curve of Fig. ||). 

Comparing incident photon energies just above and below the 33.17 keV K- 
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FIG. 2. Simulated instrumental energy 
loss spectra, commonly known as count 
spectra, for monoenergetic photons in- 
cident on a Spectroscopy Detector, the 
detector's surroundings on CGRO and 
the Earth's atmosphere. The simulation 
matches the conditions with which SD 
5 observed GRB 920311 (BATSE trig- 
ger 1473), except for clarity the simu- 
lated data extend to energies slightly be- 
low where the current calibration is reli- 
able. (Data from SD 5 and GRB 920311 
are shown in Fig. M.) In each case the 
input flux is 1 photon s~^ cm~^. Left 
curve: photon energy 20 keV; middle, 
bold curve: 50 keV; right curve: 500 keV. 



edge of iodine, the higher energy photon has a higher interaction probability, 
but more importantly, a lower probability of complete energy absorption. For 
a hard incident spectrum, such as a GRB spectrum, this K-edge effect results 
in a deficit from 33 to w 50 keV in the count spectrum, which appears as a 
peak at ~ 30 keV and is predicted by the detector response model (see Fig. 0). 
This feature has been used as a verification of the calibration of the SDs ( |3^ ) 
and should not be mistaken for an astrophysical line. 

At 500 keV, photoelectric absorption occurs only « 20% of time, making a 
contribution to the full-energy absorption peak (see Fig. |2[ right curve) . The 
most likely interaction is Compton scattering, which transfers only a portion of 
the incident photon's energy to an electron. Complete energy absorption will 
occur only if additional interactions occur, either photoelectric or Compton. 
In Compton scattering the maximum energy transfer to the electron, 331 
keV, occurs when the incident photon is scattered 180°. A range of scattering 
angles approaching 180° creates the peak in the curve at 310 keV and the lack 
of larger energy transfers in single scattering events causes the valley above 
331 keV (see Fig. ^). Correspondingly, the minimum energy of a scattered 
photon is 169 keV, so a range of angles approaching 180° for incident photons 
that scatter into the detector from the spacecraft or the Earth's atmosphere 
creates the peak at 190 keV. 

At yet higher energies pair production becomes important. 

SPECTRAL DECONVOLUTION AND LINE DETECTION 

As described in the previous section, an observed 20 keV count could be 
due to a 20, a 50, or even a 500 keV photon! When a single count is observed, 
it is impossible to deduce the energy of the incident photon. If many counts 
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are observed, the incident spectrum can be deduced within statistical limits 
in the process known as deconvolution. 

We approximate the continuous (as a function of energy) detection process 
with the following discrete equation: 

c = Bp. (1) 

Here c is the vector with the counts versus energy bins. Normally these data 
are only available in binned form in order to reduce telemetry requirements — 
as long as the energy bin widths are small compared to the detector's energy 
resolution essentially no information is lost. Similarly, p is the vector repre- 
senting the incident photon spectrum. The detector is represented by D, the 
detector response matrix, which is obtained via Monte Carlo simulations of 
gamma-ray interactions in a computer model of the detector (e.g., (|35|)). 
The obvious solution is 

p = T)-'^c. (2) 

This simple approach does not work if one wishes to deduce information at a 
resolution at or better than the intrinsic resolution of the detector, which is 
the goal of analyzing data for the presence of lines. The problem is that if the 
energy widths of the rows and columns of D are comparable to the detector 
resolution, then neighboring columns will be very similar and D will be nearly 
singular. Inverting D will be numerically unstable and the solution p will be 
unreliable, especially in the presence of statistical fluctuations in the observed 
count spectrum Cobs- 

Instead, the standard approach in astrophysics has become that of forward- 
2|). A parameterized spectral model is assumed and used to 
calculate Pmodci, from which eq. 1 yields a model count spectrum Cmodci- The 
model count spectrum Cmodei is compared to the observed count spectrum Cobs 
using some statistical measure such as or likelihood and the parameters of 
the photon model are optimized so as to minimize the discrepancy between 
Cmodei and Cobs according to the chosen statistical measure. 

It is very important to realize that a solution obtained by forward-folding 
is not unique but is rather photon model dependent ( pl] , p5| ). Even if a model 
fit results in a good value, another, possibly unknown, model might result 
in an equal or better value. This is shown in Fig. ^, where solutions 'a' 
and 'b', based upon different photon models, have essentially identical 
values. What a forward- folding solution provides is the parameter values of 
the assumed photon model ( ^5|) . 

In gamma-ray astrophysics, deconvolved photon points are frequently ob- 
tained by scaling the model photon spectrum by the ratio of the observed 
over modeled counts (p5| ) — Fig. || (right) is an example. This practice is po- 
tentially misleading because the deconvolved points are model dependent yet 
there is an almost irresistible temptation to regard them as incident spectrum 
"data" points. In the examples in Fig. ^, based upon alternative continuum 
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FIG. 3. To illustrate the photon model dependent nature of the forward-folding 
deconvolution method, identical data is analyzed using four different photon mod- 
els. Model 'a' assumes Band's spectral form obtaining •x^=22'i for 200 de- 
grees-of-freedom (dof). Model 'b' is a power-law with two breaks; x^=227 for 198 
dof. Model 'c' is a power law times an exponential cutoff; x^=244 for 201 dof. Model 
'd' is a power-law with one break; x^=267 for 200 dof. The data are for the interval 
1.088 to 26.624 s of GRB 920311 as observed by SD 5; above 200 keV, the data 
have been rebinned into wider channels for display purposes. Each group after 'a' 
is shifted downwards by xvTO. Left: count data (points) and models (histograms). 
Right: Deconvolved points and models. The uFi, spectrum is E^x the photon flux 
spectrum. The photon models obviously differ — there are also differences in the 
values of the deconvolved points, e.g., the values of right-most point are 483 ± 101, 
483 ± 101, 593 ± 113 and 449 ± 98 keV s"^ cm'^ for models 'a', 'b', 'c' and 'd', 
respectively. 



models, the differences in the deconvolved points are subtle and the practice 
not too pernicious. When a line with width comparable to or smaller than 
the detector's intrinsic resolution is considered, the deconvolved spectra can 
exaggerate the significance of a line — see Fig. |[ Papers analyzing spectra in 
the X-ray band generally show a graph of the count data and model and a 
graph of the residuals (e.g., (p2|)). Only rarely is a graph of a deconvolved 
spectrum presented. The gamma-ray community is advised to emulate this 
practice. 

Because of the limitations of gamma-ray spectral data, the following steps 
are necessary to demonstrate that a line feature exists: 

• Deconvolve spectra using the forward-folding technique. This is the stan- 
dard approach in X-ray and gamma-ray spectroscopy. 
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FIG. 4. To illustrate the dependence 
of the deconvolved photon points on the 
model for a model containing a line, 
two deconvolutions of the same data are 
shown. For the top curve, a continuum 
model is assumed, while for the bottom 
curve (shifted downward by x2.5), an 
additive Gaussian line is also included 
in the model. The data are that of 
the line candidate in GRB 940703— for 
count data and models see ref. (Mi. 
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Show the data and the model in counts rather than deconvolved photons. 
The fit is performed by comparing observed and model counts and the 
data and model should be displayed in this space. 

Show that the line is significant versus all reasonable continuum models. 
Suppose that Ax^ indicates that a line is a statistically significant im- 
provement to continuum model A — one might be tempted to say that 
the line is "proven". However, if the line insignificantly improves contin- 
uum model B and model B has comparable to that of A+linc, then 
we cannot say the line must exist because continuum model B may be 
the correct explanation. Since we do not know the correct continuum 
model of GRBs, we must try all reasonable ones (|ll|,^2|); since all mod- 
els cannot be tried, judgement is needed. In practice, this means that 
a sufficiently flexible continuum model should be used With cur- 

rent data, an appropriate choice for all bursts is Band's function (^5), 
which is usually the best-fitting model for bright bursts (e.g.. Fig. 2). 
A power law with multiple breaks is also sufficiently flexible, but has 
the disadvantage of not being smooth. An evolving or two-component 
continuum spectrum might mimic a line when a single continuum model 
is used for the analysis (^ 36). 

Provide a quantitative evaluation of the statistical significance of the 
feature. Without such an evaluation one does not know the strength of 
the evidence that a line is real rather than a fluctuation. It has been 
traditional to use the F-test (e.g., (|l|,|,|3|)), but the BATSE team has 
recently realized that Ax^ is more appropriate for this case, where one 
knows the uncertainties from Poisson statistics (|2^,|l^,0). In evaluating 
the significance of a feature, the number of trials should be considered. 
Consider possible systematic errors. Does the line centroid match that 
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of a background line? If so, is the background subtraction adequate? Is 
the detector model and calibration adequate? 
• If possible, evaluate the feature for confirmation/consistency. If more 
than one instrument, or more than one detector of an instrument, ob- 
served the burst with appropriate energy coverage and resolution, then 
the two observations should be compared. Ideally, both detectors would 
show the feature to be statistically significant, thereby confirming the 
detection. Depending on the sensitivity of the detectors and the strength 
of the burst, this may not be possible. In any case, consistency between 
the data of the two detectors is required. 

LINE OBSERVATIONS 

Table 1 summarizes the analyses methods used in previous observations of 
gamma-ray burst spectral features below 100 keV and Table 2 presents the 
results of those analyses. Many of the previous analyses were done before the 
gamma-ray community realized the importance of all of the steps listed in 
the previous section and before computer improvements made the forward- 
folding approach practical. Probably due to improvements in the sensitivity 
of detectors, there has been historical progression towards continuum models 
with more parameters (^. Many of the earlier line analyses were done using 
continuum models that are now known to be too simple, raising the possibility 
that the line detections are an artifact of the chosen continuum model. Unfor- 
tunately, these problems mean that some of the previous analyses must now 
be judged as inconclusive with regard to the existence of gamma-ray burst 
lines. Considering the importance of this question, reanalyses are merited. 

The Konus observations pioneered GRB line studies. The analysis tech- 
nique was cleverly designed to minimize the computations required: the count 
spectra were deconvolved using a standard model and the resulting photon 
spectra were iteratively improved (^8|). The continuum model used was op- 
tically thin thermal bremsstrahlung, oc exp ~E / Eq. The explanation of 
the deconvolution procedure shows one count spectrum as an example (|2^). 
It has been frequently suggested that the Konus lines might be harmonics of 
fundamentals that are below the typical detector threshold of 30 keV. 



TABLE 1. Low-Energy Spectral Features: Analysis Methods 



Instrument 




References 


Forward- 
Folding 
used? 


Count Data 
& Model 
shown? 


Flexible 
Continuum 
Model used? 


Statistical 
Significance 
evaluated? 


Konus 
HEAD A-4 
Ginga 
Lilas 


( 
( 
( 
( 


26 

H 

12 

f 


m 


No 
No 
Yes 
Yes 


No 
No 
Yes 
No 


No 
No 
Yes 
Yes 


No 
Yes 
Yes 
Yes 
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TABLE 2. Low-Energy Spectral Features: Results 



Instrument 


GRB 


References 


A 2 

Ax 


Centroids 


equiv. width 










(keV) 


CkeV") 


Konus 


many^ 


— r-rrr 




30-70 


5-15 


HEAO A-4 


780325*^ 


(|i|) 


14.0 


76±5 


24 




780325'^ 




16.6 


49±3 


13 




780608 






66±7 


14 


Ginga 


880205 




48.5 


19.7±0.7, 38.0±1.6 


3.7, 9.1 




870303=* 




30.6 


21.1±1.1 


10.5 




870303° 




26.2 


21.4±0.7, 2x21.4 


4.8, 8.4 




890929 




25.6 


26.3±1.5, 46.6±1.7 


4, 8 


Lilas 


890306 




89.0 


11.2±0.5, 34.6±1.4 





^Unique amongst the results listed in this table, one GRB had an emission feature. 
''First peak of the burst. 
'^Second peak of the burst. 
•'Interval 'SI'. 
•^Interval 'S2'. 



The most significant line in the HEAO A-4 data is in the 2nd peak of 
GRB 780325 (|9|). The line significance was evaluated using a simple contin- 
uum model, an exponential, fit to only the data below 200 keV ( |l9| ) . Using a 
better model and all of the data might raise or lower the significance. The de- 
convolution was done using photons-to-counts efficiencies derived for an E^'^ 
spectrum, a procedure equivalent to approximating D with a diagonal matrix. 
This procedure was justified by the suppression of the Compton scattering re- 
sponse by the active shielding of the instrument and has the advantage of 
conservatively assuming that lines do not exist. The feature is stated to be 
visible in two detectors, but only summed data are shown (p^. 

The analysis of the Ginga data is excellent. (However, some papers present 
only deconvolved spectra.) The chance probability of obtaining a reduction 
in of 48.5 by adding two lines, assuming that the lines do not exist, is 
10"^! The harmonic relation between the line centroids is powerful evidence 
for the cyclotron resonant scattering interpretation. A possible concern is that 
the lines were in the same channels in all 3 bursts (^3|) — while the centroid 
energies are different in GRB 890929, for this event the gains of the detectors 
were lower than normal (^). 

The recent results of Lilas on GRB 890306 are even more statistically 
significant — lines at 11 and 35 keV are reported in a 68 s interval with an 
F-test probability of 2 x lO"'^^ that this is a chance result (H). The value 
for the model containing the lines is somewhat high: 43.2 for 26 degrees-of- 
freedom. This value, along with the location of one line near the detector 
threshold and the other near the iodine K-edge, are areas for concern. It 
would be desirable to see the count model for these data and also for a com- 
parable spectrum not containing lines. The spectral evolution during the 68 s 
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interval, which is essentiaUy the entire burst, should also be investigated. 

BATSE AND LINES 

The Spectroscopy Detectors were added to BATSE because of the Konus 
line observations. The pre-launch expectations of the BATSE team were that 
lines would be readily found. Simulations and performance tests ( |3^ ) 
verify the ability of BATSE to detect Gm^a-like lines, with generally poorer 
sensitivity than Ginga for 20 keV lines and better for 40 keV lines. However, 
no lines have been detected with BATSE (^,^. While disappointing, this 
lack of line detections by BATSE is not yet in strong contradiction with the 
observations of Ginga (pSPJ^). 

BATSE has a number of advantages for line studies: 

• State of the art scintillation detectors: excellent energy resolution and 
advanced electronics incorporating active PMT bleeder strings and base- 
line restoration to handle large pulses and high count rates, 

• Excellent temporal resolution so that lines can be found on whatever 
time scale they may exist, 

• GRB locations from the LADs to aid analysis of the SD data, 

• Good sensitivity, 

• Extensive performance verification, on the ground and in-orbit (^2[), 

• Multiple detectors and detector types, which enable consis- 
tency/confirmation studies. 

The BATSE team is therefore confident of the instrument's ability to de- 
tect lines. Because of the lack of detections by the initial approach, visually 
examining burst spectra for lines (p3|), the team has implemented the most 
thorough line search ever conducted. This comprehensive computer search 
examines essentially all time scales and energy centroids below 100 keV. The 
new approach has already identified 8 candidate features with Ax^ > 20 (||) . 
Currently analysis is in progress to determine the consistency of the multi- 
detector data and hopefully to confirm some of the features with multiple 
detectors. Until that work is completed the BATSE team considers these 
features to be candidates rather than detections. 

THE FUTURE 

There is much more work to be done. On the interpretational side, bursts 
from galactic halo or cosmological distances have intrinsic luminosities 10^ or 
more greater than deduced in the days of the galactic disk paradigm. Con- 
sequently, the physics of line creation will be different than previously en- 
visioned. This has been examined in the context of sources in the galactic 
halo (e.g., (|2|)), but little work has been done in the context of cosmological 
models. Several of the past line observations merit reanalysis using techniques 
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developed since they were published. The analysis of the BATSE data con- 
tinues in order to determine the reality of the candidates identified by the 
comprehensive search. With eight candidates already identified, and typically 
several high-gain SDs viewing a burst, the BATSE team will be able to deter- 
mine the consistency of the data and confirm many of the candidate features 
or to demonstrate that there is some problem. Data are also currently being 
collected and analyzed with the scintillation detectors of the Konus-W instru- 
ment (0,^9|) and the high energy-resolution detectors of TGRS ( |3^ ) . The near 
future will bring the launch of HETE ^ and PGS on Mars-96 (|o|). 
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